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We present a determination of the Hubble constant from measurements of the 
Sunyaev-Zeldovich Effect (SZE) in an orientation-unbiased sample of 7 2 < 0.1 
galaxy clusters. With improved X-ray models and a more accurate 32-GHz cal- 
ibration, we obtain Ho = 64^ J* ± 14 sys kmsec — 1 Mpc - 1 for a standard CDM 

cosmology, or Ho = 66^** ± 15 S yskmsec _1 Mpc" 1 for a fiat ACDM cosmology. 
The random uncertainty in this result is dominated by intrinsic CMB anisotropics, 
and the systematic uncertainty is chiefly due to uncertainties in the X-ray spectral 
models. The implications of this result and future prospects for the method are 
discussed. 



1 Introduction 

For over two decades it has been known that the combination of X-ray and SZE 
observations of rich galaxy clusters, under the assumption of spherical symmetry, 
yields a direct measurement of the cosmic distance scalcEra. This comes about due 
to the fact that the X-ray brightness of a cluster depends upon a different power of 
h than the SZE (which is unaffected by cosmological surface brightness dimming), 
allowing the radio decrement, AT D b s = q, to be predicted from the X-ray data up 
to a power of h: AT pre d = phr 1 / 2 . Then h is given by 




Only in the past few. Kfars have reliable and accurate applications of this method 
become possible oa&Elll Due to the assumption of spherical symmetry in this 
analysis, selection biases have been a great concern in using the SZE as a distance 
measure. To address this concern Myers et al. (1997) defined an X-ray flux-limited 
sample of 11 z < 0.1 clusters and began a campaign to measure distances to these 
clusters. With this strategy, random departures from spherical symmetry in in- 
dividual clusters will not bias the sample average for Ho. With observations of 4 
clusters Myers et al. obtain a Hubble Constant of 54 ± 14 km sec -1 Mpc -1 . The 
accuracy of this result is limited by a 7% radio calibration uncertainty and esti- 
mated 15-30% X-ray model uncertainties for each cluster. In this proceeding, we 
present more recent observations in the campaign initiated by Myers et al., and 
discuss other improvements in these results. 

Our efforts are concentrated on reducing and understanding the uncertainties in 
the X-ray models, since these have always been a limiting consideration for the SZE. 
This is facilitated by the advent of abundant X-ray data from ROSAT and other 
satellites. In addition, we have conducted an absolute radio calibration program at 
the Owens Valley Radio Observatory (OVRO) which has reduced our calibration 
uncertainties by over a factor of two. We have also quantitatively evaluated the 
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effect of CMB anisotropics on our observations. 

Throughout this proceeding we we use Ho = 100 /i _1 kmsec _1 Mpc -1 . We will 
consider two cosmologies: SCDM (ft m = 1,0a = 0), and ACDM (f2 TO = 0.3, SIa = 
0.7). 

2 Cluster Sample 

In defining a sample for SZE-based Ho determinations the most important con- 
sideration is that the sample be free of any orientation bias, such as would be 
introduced, e. g. , , by a surface brightness selection. This is especially important 
due to the fact (Eq ^) that systematic errors in the observations enter quadrati- 
cally into the final Ho result. Since X-ray cluster catalogs are most reliable at high 
flux levels — and the optical catalogs which underlie most current X-ray cluster 
catalogs are most reliable at low redshifts — a low-z sample is indicated. This has 
the additional advantage that the X-ray models will be better since, at low-z, the 
signal is stronger and angular resolution is less of an issue. 

Mason & Myers (2000) define an expanded, X-ray flux-limited cluster sample of 
z < 0.1 .objects. This is a 90% volume-complete sample selected from the XBAC 
catalog i3 with the restrictions Fx > 1.0 x 10~ n erg cm~ 2 s -1 . At z = 0.1, the 
volume-completeness criterion corresponds to Lx > 1.13 x 10 44 ft.~ 2 ergs -1 (0.1 - 2.4 
keV). The resulting set of 31 clusters contains the Myers et al. sample. Although the 
XBAC sample is ultimately derived from an optical catalog, the bright, rich clusters 
which pass our luminosity criterion should be reliably selected in any orientation. 
The 7 clusters for which we present measurements here are all members of the 
smaller Myers et al. sample, but this work is part of an ongoing project to survey 
distances to the objects in our expanded sample. 

3 X-ray Models 

Mason & Myers (2000) also present X-ray models for the 22 of the 31 clusters which 
had public ROSAT PSPC data as of May 1999. The primary focus of this analysis 
was to quantify the uncertainty in the beam-convolved inverse Compton optical 
depth, t sw . To do this, the ROSAT PSPC data were smoothed, and 10 3 realiza- 
tions of each cluster were generated including Poisson noise; each realization was 
analyzed just as we analyzed the real data. We used a similar strategy to quantify 
errors induced by inaccurate vignetting correction and background subtraction. We 
find, overall, that the beam-convolved inverse-Compton optical depths are robustly 
constrained by the ROSAT data, and unaffected by realistic PSPC systematics. As 
an example, Figure |l| shows the 10 3 model predictions for the beam-convolved op- 
tical depths (t sw ) from our analysis of the A2142 data. This is an extreme example 
since we have excised the central region of the PSPC data, which is contaminated 
by cooling flow emission, and allowed the fit to be completely unconstrained there. 
The result is that the core radius o is poorly determined; nevertheless the overall 
optical depth is well constrained since the SZ decrement is mostly determined by 
the extended lines of sight. The cancellation of the large uncertainties in the in- 
dividual parameters ((3,0 o ,n eo ) comes about due to strong parameter correlations, 
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Figure 1. Predictions of the 5.5-meter beam-convolved SZE decrement for A2142. The solid line 
shows the distribution mean and the dashed lines show two-sided 68% confidence regions. 



Table f . Cluster Model Parameters 



Cluster z 




(') 



(10~ 3 h 1 ' 2 cm- 3 ) (kcV) 



A399 

A401 

A478 

A1651 

Coma 

A2142 

A2256 



0.0715 
0.0748 
0.0900 
0.0825 
0.0232 
0.0899 
0.0601 



4.33 ±0.45 
2.26 ±0.41 
1.00 ±0.15 
2.16 ±0.36 
9.32 ±0.10 
1.60 ±0.12 
5.49 ±0.21 



0.742 ± 0.042 
0.636 ±0.047 
0.638 ±0.014 
0.712 ±0.036 
0.670 ± 0.003 
0.635 ±0.012 
0.847 ±0.024 



3.23 ±0.18 7.0 ±0.2 

7.90 ±0.81 8.0 ±0.2 

27.81 ±9.7 8.4 ±0.7 

7.14 ± 3.20 6.1 ±0.2 

4.52 ±0.04 9.1 ±0.4 

14.95 ±1.0 9.7 ±0.8 

4.08 ±0.08 6.6 ±0.2 



and it is generally important to account for these {e.g., by Monte Carlo methods) 
in order to understand the uncertainties in the observables. 

A summary of the models for the 7 clusters on which we have radio observa- 
tions is presented in Table |l|. For Coma, following Hughes, Gorenstein, & Fab- 
ricant (1988), we adopt a hybrid temperature profile which is isothermal inside 
of 500/i _1 kpc and adiabatic (7 = 1.5) outside of this. For the other clusters we 
assume an isothermal profile. We use the global, cooling-flow corrected electron 
temperatures of Markevitch et al. (1998). Our value for Ho is not sensitive to 
the temperature profile we assume: it changes by < 5% if we assume hybrid pro- 
files for all of the clusters. While more extreme temperature profile models would 
have a stronger affect on our result, such models are not motivated by current 
analysesOE 3 ! 

4 Radio Observations 
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Table 2. OVRO SZE Observations 



Cluster y obs 

(io- 5 ) 




A399 3.24 ±0.41 

A401 6.93 ± 0.46 

A1651 4.88 ±0.59 

Coma 6.38 ±1.01 

A2142 9.10 ±0.52 

A2256 5.00 ±0.60 



2.5 
5.2 
5.6 
4.0 
6.7 
3.6 



4-1 Recalibration 

In March through May of 1998 we undertook a program of absolute calibration 
with the OVRO 1.5-meter radio telescope. The gain of this antenna was measured 
by NIST's Boulder, CO near-field test range in May of 1990; these measurements 
in conjunction with accurate total power determinations on-site at OVRO allowed 
a direct determination of the flux densities of celestial sources. The primary target 
for this program was the supernova remnant Cassiopeia A (Cas A), the second 
brightest extrasolar source in the sky at cm wavelengths. We find an epoch 1998 
flux density ScasA,WS8 = 194 ± 4 Jy at 32.0 GHz. Observations of Jupiter relative 
to Cas A with the OVRO 5.5-meter telescope give T Jup = 152 ± 5K at 32.0 GHz. 
All of our SZ measurements are referred to this flux density scale, which improves 
by a factor of two over the uncertainty present in the Myers et al. calibration. 
The absolute calibration measurements are discussed at length in Mason et al. 
(1999), and..are the basis of the flux density .scales which have, been adopted hx 
SaskatoonBO, the 30 GHz channel of MATtJ, and the CBI0; the RING5M0 
calibration is comparable but independent. 

4-2 SZE Observations & Effect of the Intrinsic Anisotropy 

From October 1996 through March 1998 we observed A399, A401, and A1651 with 
the OVRO 5.5-meter telescope. The results of these_pbservations, together with 
measurements of ComaQ and A2142, A2256, and A478B, are shown in Table §. The 
errors reported in this table are the observational noise only. The signficance of the 
detection a is also reported, including the effects of intrinsic CMB anisotropies. 

Intrinsic anisotropies on this scale [I = 590 ± 200) are known to be significant, 
although these will be averaged down somewhat by our scan pattern. To evaluate 
the impact of these anisotropies, we generated 10 3 realizations o£4deg 2 patches of 
sky from a ACDM power spectrum normalized to the RING5ME3 observed power 
level. Each realization was convolved with the 5-m main beam and the switching 
pattern characteristic of a typical SZE observations. We find a residual CMB 
signal of cr cm b = 64 /iK, or cr cm b = 1.24 x 10~ 5 in terms of the Compton-y, which 
is accounted for in our Hq analysis. 
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5 Interpretation &; Discussion 

When the effects of intrinsic CMB anisotropy are included, faint clusters like A399 
are seen to be detected at less than 3a. It is important to keep such clusters in the 
analysis, however, in order that an orientation bias not be introduced into our result. 
To account for these low signal-to-noise measurements, we have developed a maxi- 
mum likelihood method of estimating Ho from the SZ and X-ray data. Preliminary 
results from the application of this method to our 7 clusters are shown in Table [| 
The calibration uncertainties are 3% (radio) and 8% (X-ray), and we estimate a 
10% uncertainty in the SZE predictions due to the possibility of substructure and 
non-isothermality in the ICM. We find H = 64l^±14 SJ/s kmsec- 1 Mpc" 1 (SCDM), 
or Ho =j— £36^^ ± 15 sys kmsec _1 Mpc -1 (ACDM). Final results on this sample are 
in pressed. 

The ages implied [(10.2 ± 3.2) x 10 9 years SCDM; (14.2 ± 4.5) x 10 9 years 
ACDM] are consistent with recent age determinations from main-sequence fitting, 
which give ages of (> 12 ± 1) >j-|10 9 years Eircj. From a hydrostatic analysis of 
22 X-ray clusters in our samples, together with this value for H , we determine 
£Im = 0.33 ± 0.05. This argugs. strongly against SCDM cosmologies, and in light 
of recent CMB measurementsEZlliSEj which indicate f2 tot ~ 1, may suggesi_.another 
form of energy density such as the cosmological constant or quintessence! 2 -! These 
results are broadly in agreement with the concordance modelscj suggested by a 
number of independent means. 

CoAffirage of clusters in this sample is continuing with the Cosmic Background 
ImagerEjO (CBI). The CBI's greater coverage (~ 20 clusters) will significantly 
reduce uncertainties due to intrinsic CMB anisotropics. Although the results we 
report here incorporate improved X-ray models, the modelling uncertainties will 
still be the limiting consideration for larger samples. In order to achieve a 10% 
determination of Ho from the SZE, a better understanding of the spectral models 
for the cluster atmospheres will then be essential. 



Table 3. Ho results on 5 Clusters 



Cluster 


y bs 


Upred 


Ho 




(io- 5 ) 




(km sec -1 Mpc -1 ) 


A399 


3.24 ±0.41 


3.27 ±0.20 


102±^ 


A401 


6.93 ±0.46 


4.82 ±0.32 


48±?| 


A478 


7.77 ±0.58 


6.05 ±0.54 


fil +33 
D1 -20 


A1651 


4.88 ±0.59 


2.92 ±0.17 


36±?| 


Coma 


6.38 ± 1.01 


5.00 ±0.38 


62±g 


A2142 


9.10 ±0.52 


8.12 ±0.74 


70+34 

<y_24 


A2256 


5.00 ±0.60 


4.11 ±0.26 


D( -28 


SAMPLE 









NOTE: Observational and model ncertainties are la random errors only. 
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